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Abstract 

We consider fermions on an extra dimensional interval. We find the boundary condi- 
tions at the ends of the interval that are consistent with the variational principle, and 
explain which ones arise in various physical circumstances. We apply these results 
to higgsless models of electroweak symmetry breaking, where electroweak symmetry 
is not broken by a scalar vacuum expectation value, but rather by the boundary 
conditions of the gauge fields. We show that it is possible to find a set of boundary 
conditions for bulk fermions that would give a realistic fermion mass spectrum with- 
out the presence of a Higgs scalar, and present some sample fermion mass spectra for 
the standard model quarks and leptons as well as their resonances. 



1 Introduction 



The most exciting question facing particle physics is how electroweak symmetry is broken 
in nature. Since the scattering of massive W and Z bosons violate unitarity at the scale of 
~ 1.8 TeV, we know that some new particles must appear before those scales are reached to 
unitarize these amplitudes or the theory will be strongly interacting. In 4D the only possi- 
bility to unitarize these scattering amplitudes (with a single particle) is via the exchange of 
a scalar Higgs particle. It has been recently pointed out in jTJ that extra dimensions may 
provide an alternative way for unitarizing the scattering of the massive gauge bosons via 
the exchange of a tower of massive Kaluza-Klein (KK) gauge bosons (see also in this 

case electroweak symmetry would be broken not by the expectation value of a scalar (or 
a scalar condensate), but rather by the boundary conditions (BC's) for the gauge fields.* 
As long as the BC's are consistent with the variation of a fully gauge invariant action, 
the symmetry breaking will be soft in the sense that the UV properties of the scattering 
amplitudes will be as in the higher dimensional gauge theory pQ. 

A model of higgsless electroweak symmetry breaking (EWSB) with a realistic gauge 
structure has been presented in jB]. However, the Higgs scalar of the standard model (SM) 
serves two purposes: besides breaking the electroweak symmetry it is also necessary for the 
generation of fermion masses without explicitly breaking gauge invariance. The purpose of 
this paper is to examine how fermion masses can be generated in higgsless models where 
EWSB happens via BC's in extra dimensions. 

The structure of these higgsless models is generically of the following form :6J: we 
consider the a modification of the Randall- Sundrum model jH] with gauge fields in the 
bulk PHI]], where the bulk gauge group is SU(2) i xSU(2) R xU(l) J B_z / . The addition of the 
second SU(2)# in the bulk is necessary [TT] in order to ensure the presence of a cus- 
todial SU(2)^ symmetry in the holographic interpretation [T2|. On the Planck brane 
SU(2)^x U(1)b-l is broken to U(l)y, while EWSB happens on the TeV brane where 
SU(2) L xSU(2) fi — >SXJ(2) D (an early flat space version of this model was presented in [I], 
and recently re-examined in ^3]). Two important issues regarding the warped higgsless 
model that were not addressed fully in [U] were the generation of fermion masses without a 
Higgs, and the corrections to electroweak precision observables. The issue of fermion masses 
will be discussed in detail in this paper. There are several potential sources for corrections 
to electroweak precision observables in a higgsless model: the enlarged gauge structure, 
the missing Higgs scalar, and the modified fermion sector. Recently ^3] examined the S 
parameter in the flat space version of this model and found that analogously to technicolor 
theories there is a large positive contribution. However, no comprehensive analysis of the 
electroweak observables including all sources of corrections listed above has been done to 
date. We plan to address these issues for the case of an AdSs bulk (as considered in jHllZl) 
in a forthcoming publication. 

*For other possibilities of utilizing extra dimensions for electroweak symmetry breaking see 
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In order to be able to generate a viable spectrum and coupling for the SM fermions, 
the fermions have to feel the effect of EWSB, so they need to be connected with the TeV 
brane. However they can't be simply put on the TeV brane, since in that case they would 
form multiplets of SU(2)£>, which they don't. Thus the fermions also have to be put into 
the bulk, as in |14[ll5|ll6|ll7|ll8|ll9j . We assume that the left handed SM fermions will form 
SU(2)l doublets and the right handed ones SU(2)# doublets (including a right handed 
neutrino). Since 5D bulk fermions contain two 4D Weyl spinors (like a 4D Dirac fermion), 
one has to first make sure that in every 5D bulk fermion there is only a single 4D Weyl 
spinor zero mode. These zero modes will be identified with the usual SM fermions. In 
order to recover the usual gauge coupling structure for the light fermions, the zero modes 
for the light fermions have to be localized close to the Planck brane. Since the theory 
on the TeV brane is vector-like, one can simply add a mass term on the TeV brane that 
connects the left and right handed fermions. This is however not sufficient, since the gauge 
group on the TeV brane would force the up-type and down-type fermions to be degenerate. 
The splitting between these fermions can be achieved by mixing the right handed fermions 
with fermions localized on the Planck brane where SU(2)# is broken (which is equivalent 
to adding different Planck-brane induced kinetic terms for the right handed fermions). 
The detailed models for the fermion masses for the warped space higgsless model will be 
presented in Section 7 (while the analog constructions for the somewhat simpler flat-space 
toy model can be found in Section 5). 

Before we discuss fermion mass generation for the higgsless models in detail, we will 
discuss the general issues surrounding the often confusing subject of BC's and masses for 
fermions in one extra dimension. In Section 2 we examine the possible BC's for fermions on 
an interval that are consistent with the vanishing of the boundary variations of the action. 
This is an extension of the general discussion of |I] of BC's in an extra dimension to the 
fermion sector. In Section 3 we discuss the general KK decomposition for fermions on an 
interval (with some simple examples worked out in detail in Appendix B), while in Section 
4 we give the physical interpretation of the various BC's obtained from the variational 
principle. More important examples for BC's in the presence of mixing of bulk fermions on 
a brane are presented in Appendix C. In Section 5 we apply the results of Sections 2-4 to 
propose the BC's for the fermion sector of the flat space higgsless model. In Section 6 we 
discuss the general issues of fermionic BC's in warped space, and then finally present the 
BC's and mass spectra for the warped space higgsless model in Section 7. We conclude in 
Section 8. Appendix A contains notations and spinor conventions. 



2 Fermion Boundary Conditions from the Variational 
Principle 

We start by considering a theory of 5D fermions on an interval of length L, with a bulk 
Dirac mass m, and possibly also masses for the component fermions on the boundaries. 
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For the moment we will assume that the geometry of the interval is flat and we will come 
back later to the phenomenologically more interesting case of an AdSs interval. In 5D, the 
smallest irreducible representation of the Lorentz group is the Dirac spinor, which of course 
contains two two-component spinors from the 4D point of view. The bulk action for the 
Dirac spinor \& is given by the usual form 

S = j d 5 x Y M d M ^ - <3m* T M ^) - mtftf J (2.1) 

where M = 0,l,2,3,5. Usually, in the second term the differential operator is integrated by 
parts and gives a contribution identical to the first term. However when the fifth dimension 
is a finite interval, boundary terms appear in the process of integrating by parts and using 
the conventional form of the action would require us to explicitly introduce those boundary 
terms. That this is the convenient starting point can be seen from the fact that the action 
is real or equivalently that the corresponding Hamiltonian is hermitian. Writing out the 
action in terms of the two-component spinors contained in the 5D Dirac spinor as 



* = £ ( 2 - 2 ) 




and integrating by parts in the 4D coordinates where we do require that the fields vanish 
at large distances, we obtain the following Lagrangian for the two-component spinors (for 
spinor and gamma matrices conventions, see Appendix EJ: 

S= I d 5 x {-ix^d^x ~ iip^d^ + \ {^% X ~ X%$) + m(ij>x + xVO) , (2.3) 



where d$ =85 — 85, with the arrows indicating the direction of action of the differential 
operator. Varying the action with respect to x an d ip we obtain the standard bulk equations 
of motion which are given by 

— io^d^x ~ ds' 1 !' + m i' — 0, 
-io^d^i) + d 5 x + mx = 0. (2.4) 

However one needs to be careful with the variation, since one needs to do an integration by 
parts in the extra dimension, which will give an extra term for the variation of the action on 
the boundaries of the interval. Requiring that the boundary term in the variation vanishes 
will give the desired boundary conditions for the fermion fields (we denote by [X]q the 
quantity X\ L — X ): 

SS bmlnd =~ J d 4 x [5x^P - Sifix ~ S$X + <5x4>]o = 0. (2.5) 

To proceed further we need to specify the boundary conditions. These have to be 
such that the boundary variation in (|2.5J) vanishes. Note, that this is a somewhat unusual 
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boundary variation term (at least compared to the case of scalar and gauge fields) since it 
mixes the two Weyl spinors. We will first discuss the simplest and most commonly adopted 
solutions, and then consider the more general cases. The most obvious solution to enforce 
the vanishing of (|2.5jl is by fixing one of the two spinors to zero on the endpoints, for 
example 

V>|o,z = 0. (2.6) 

As soon as we imposed this condition, we also have that 5i/j\o,l = 0, and the full boundary 
variation term vanishes. This would naively suggest that \ remains arbitrary at the end- 
points, however this is not the case, since we still have to require that the bulk equations of 
motion are satisfied everywhere, including at the endpoints of the interval. Since the bulk 
equations mix i[) and x, when ip = we get a first order equation for just x> which can be 
considered as the boundary condition for the \ field: 

(d 5X + m X ) l0iL = 0. (2.7) 

In the limit of m — > this is the BC that is usually employed when considering orbifold 
projections. The usual argument is that if one assigns a definite parity to x an d under 
y — > —y, then due to the term ipd$x i n the bulk, x an d have to have opposite parities, so 
if ip is chosen to have negative parity (that is it vanishes on the endpoints) then x has to 
be positive, so its derivative should vanish. This is basically what we see in this simplest 
solution, except that it is very easy to deal with the bulk mass term. If one were to think 
of this interval as the orbifold projection of a circle, then the only way one can fit a bulk 
mass into the picture is if the bulk mass is assumed to switch signs at the orbifold fixed 
points (the mass itself has negative parity), which then makes figuring out the right BC's 
in the presence of the mass term quite hard. We can see that in the interval formulation 
one does not have to worry about such subtleties. 

We have seen above that the simplest possible solutions to the vanishing of the bound- 
ary variation leads to the boundary conditions generically employed when considering orb- 
ifold constructions. However, one does not need to require the individual terms in ()2.5|) 
to vanish, it is sufficient for the whole sum to vanish. In fact, requiring the individual 
variations to vanish over-constrains the system, as is clear from a simple counting of the 
degrees of freedom of the theory. There are two constants associated with the solutions of 
two first order differential equations. One boundary condition at each end of the interval 
then specifies the system. If one forces the individual terms in the boundary variation on 
one endpoint to vanish, then there is no freedom of boundary conditions on the opposite 
endpoint. Thus one should generically only impose one BC at each end of the interval. 
Such a BC expresses one of the spinors in terms of the other. With that in mind we can 
see that the most general solution to the vanishing of the boundary variation is when, on 
the boundary, the two fields ip and x are proportional to each other: 

^|o,l= (Mjxfi + N^) (2.8) 
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where M and N are two matrices that may involve some derivatives along the dimensions 
of the boundary. The action will then have a vanishing boundary variation provided that 
M and N satisfy the two conditions 

Mj = a M e^M^e Sa and N a ^ = a N N^ (2.9) 

where ctm,n are signs due to some possible integration by part of the differential operators 
contained in M and N. Two simple solutions are 

Mj 3 = c5oP for any constant c, (2-10) 
N a p — i° ®n for any real constant c. (2-11) 

Note that for fermions belonging to a complex representation of the gauge group, gauge 
invariance requires that either the operator M or its inverse vanishes. Let us discuss in 
more detail the solutions of the type ip\o,L = Co,lX\o,l, for arbitrary values of Co,£. A 
better way of expressing this condition is by saying that some linear combination of the 
two fermion fields has to satisfy a Dirichlet boundary condition on both ends. However, 
these can be different combinations on the two sides: 

so,l i/>\o,l + Co,L X\o,l = (2.12) 

where so,l (co,l) stand for the sine (cosine) of some (possibly complex) angles, «o,l, that 
determine which linear combination of the fields on the two boundaries are vanishing. If 
there are gauge symmetries in the bulk, under which the fermions transform, then ip and x 
transform in complex conjugate representations, as can be seen from Eq. (|2.2|) . This means, 
that it is only possible to mix the two fields on the boundary with non-trivial angles if the 
fermion is in a real representation. Thus for real representations sq,l could in principle be 
arbitrary, however for complex representations the only possibilities are so,l = or so,l = 1- 
We will see later on that this choice of BC's corresponds to adding a Majorana mass on 
the brane. 



3 Kaluza-Klein Decomposition 

Now we would like to discuss how to perform the Kaluza-Klein decomposition of these 
fields. In general, when the fermion belongs to a complex representation of the symmetry 
group, the KK modes can only acquire Dirac masses and the KK decomposition is of the 
form 

x = ^2g n (y)xn(x), (3.1) 

n 

^ = J2fn(y)Mx), (3.2) 

n 



5 



where \ n and ijj n are 4D two-component spinors which form a Dirac spinor of mass m n and 
satisfy the 4D Dirac equation: 

- ia^d^Xn + m n 4> n = 0, (3.3) 
-ia^d^n + m n Xn = 0. (3.4) 

Plugging this expansion into the bulk equations we get the following set of coupled first 
order differential equations for the wave functions f n and g n : 

g' n + mg n - m n f n = 0, (3.5) 
fn-mf n + m n g n = 0. (3.6) 

Combining the two equations we get as usual decoupled second order equations: 

g" + [ml - m 2 )g = 0, (3.7) 
f" + (ml - m 2 )f = 0. (3.8) 

Depending on the sign of ml — m 2 the wave functions g n and f n will be either sines and 
cosines or sinhes and coshes (we define ccos k n L = cosh k n L for k 2 n = m 2 — ml > and 
ccos k n L = cos k n L for k\ = ml — m 2 > and similarly for ssin k n L): 
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g n {y) = A n ccos k n y + B n ssin k n y, (3.9) 
fn{y) = C n ccos k n y + D n ssin k n y. (3.10) 

These wave functions are analogous to the ones obtained for bosonic fields. For fermions, 
the bulk equations are going to teach us something more about the wave functions. Indeed 
the first order coupled differential equations (|3.5j) - (|3.6|) relate the coefficients A n , B n , C n , D n 
to each other. Using the form (j3.9|) - ()3.10|) of the wave functions and for m n ^ 0, the two 
bulk equations are equivalent to one another and impose, for kl = ml — m 2 > 0, 

mC n k n D n m n A n = 0, (3.11) 
k n C n + mD n - m n B n = 0. (3.12) 

When m 2 — m 2 n > 0, the sign of the term involving k n in the second equation is flipped. 

The boundary conditions may also allow the presence of a zero mode which can have a 
non trivial profile of the form (|3.9|) - (|3.1U|) with kl = m 2 for a non vanishing bulk mass. For 
the case of the zero mode the bulk equations (j3.5j) - (j3.5|) are decoupled and simply reduce 
to 

A = -Bo and C = D . (3.13) 

Some explicit examples of KK decomposition are given in Appendix IB1 when BC's of 
the form sif}\ +c\\ are imposed at + and L~ . We also discuss there how to amend the form 
of the decomposition (J3.1)) when the gauge quantum numbers of the fermion allow the KK 
modes to have Majorana masses. 
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Before we close this section, we just remind the reader what the status of all of these 
various boundary conditions with respect to each other is: that is can we have the different 
modes corresponding to different boundary conditions present in the theory at the same 
time? The answer is no. A theory is obtained by fixing the BC's for the fields (picking 
one of the possibilities from the list given above) once and for all. If we were to include 
modes corresponding to different BC's into the theory, we would loose hermiticity of the 
Hamiltonian, that is the theory would no longer be unitary. If one were to insist on 
putting these different modes together, an additional super-selection rule would have to be 
added that would make these different modes orthogonal, signaling that they belong to a 
different sector of the Hilbert space, which practically means that a new quantum number 
corresponding to the choice of the BC would have to be added. But usually this is avoided 
by simply considering only a fixed BC. 

4 Physical Interpretation of the Boundary Conditions 

We would like to have an intuitive physical picture of the various fermionic boundary 
conditions. Unlike a scalar field, which in the absence of boundary interactions naturally 
has a flat profile (<9s0 = 0) on the boundary, the fermions cannot have a purely flat wave 
function. This is a result of the dynamics of 5D fermions, which can be broken up into 
two two-component spinors. The bulk equations of motion for these spinors imply that, in 
the absence of a bulk Dirac mass (as we will assume throughout this section), if one spinor 
profile has zero derivative at a boundary, then the opposite spinor must obey Dirichlet 
boundary conditions, and vice versa. However, we have seen in Section 2 that there is 
a variety of BC's that one can impose instead of the simplest d$x\ — ^\ — condition. 
The purpose of this section (and its continuation in Appendix C) is to understand what 
physical situations the various BC's (some of which may seem quite obscure at first sight) 
correspond to. 

What we would like is to be able to consider a setup with arbitrary localized mas- 
ses/mixings/kinetic terms and translate these into some BC's similar to the form of (|2.8|) . 
However, it is not easy to arrive at these BC's from the variational principle if the localized 
terms are directly added at the boundary. The reason is that due to the first order nature 
of the bulk equations of motion the presence of a localized term necessarily implies a 
discontinuity in some of the wave functions. If the localized term is added directly at the 
boundary, one would have to treat the values and variations of the fields at the boundaries 
as independent from the bulk values which makes the procedure very hard to complete. 
Instead, our general approach to treating the localized terms will be the following: 

• Push the localized terms at a distance e away from the boundary, which implies the 
presence of a 5-function in the bulk equations of motion; 

• Impose the simplest BC's d 5 x\ — ^\ — at the real boundary y = 0, L; 
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• By combining the jump equation at y = e with the BC's at y = obtain a relation 
between the fields at y = e; 

• Take the limit e — > and treat the relation among the fields at y = e as the BC's for 
the theory on an interval with the localized terms added on the boundaries. 

By construction, the BC's obtained this way will always satisfy the variational principle 
(that is make (|2.5jl or its analog in the presence of more fields vanish), but this way the 
physical interpretation of the possible parameters appearing in the BC's will become clear. 

In this section we will first show what a possible physical realization of the usually 
applied simple BC d^x\ — ip\ — is, and then consider adding a Majorana mass on the 
boundary as discussed above. Many more important examples of adding localized terms 
will be discussed in Appendix C. 



4.1 Physical Realization of the Simplest Dirichlet— Neumann BC's 

The easiest way to realize the simple d§x\ = ip\ = boundary condition is to give the 
fermions a mass that is a function of an extra dimensional coordinate, where this extra 
dimension is infinite in extent. For example, one can construct a square well mass term 
given by 

m(y) = m_ 9(—y) + m + 9(y — L) (4-1) 

where 9 is the usual Heaviside function. Taking m_ and m + to be large, yet finite, the so- 
lution consists of modes which tail off exponentially outside of the well. The wavefunctions 
of x an d ip in region I (y < 0) and II (y > L) are given by 

{f = C~ e k " y i f = C + e k ™ y 

* = = (4 ' 2) 

where k^ 2 = m± — As m grows, the exponentials drop off more and more quickly, 

and thus the fermions are confined to a "fat brane" of width L (which corresponds to the 
well). The wavefunctions at y — and y = L are continuously matched with the solution 
within the "well" 

(0 < y < L) { fn = C A n C ° S mnV + ?T ^ (4-3) 
~ [ g n = A n cosm n y + B n smm n y. 

The existence of a 4D massless mode depends on the details of the mass profile m(y). 
For instance, when m + = m_ = m, a quick calculation shows that this particular mass 
background does not lead to a normalizable chiral zero mode profile. Indeed the "bulk" 
equations of motion together with the continuity conditions at y = and y = L leads to 
the quantization equation 

m n tan m n L = a/ m 2 — m 2 n , (4.4) 



8 



which obviously does not allow a massless mode. 

However, if one changes the mass profile to m + = — m_ = m, the quantization equation 

becomes 

A/m 2 — m 2 tan m n L = —m n , (4.5) 

which now supports a massless solution. This is a stepwise analogue of the well known 
domain wall localization of chiral fermions When m — > oo, the solution within the 
well can be equivalently obtained by ignoring the exterior regions and by imposing the 
following boundary conditions (for m > 0) 

<9 5 X|o = 0, V|o = 0, d 5X \L = 0, V|l = 0. (4.6) 



4.2 BC's in the presence of a brane localized Majorana mass 

We now consider adding localized terms to the fermion action, and ask how the simple 
BC derived above will change in the presence of these terms. We will illustrate in detail 
how to implement the steps outlined at the beginning of this section for the case when a 
Majorana mass is added on the boundary. This example in the context of orbifolds has 
been discussed in [22] (see also [2H])- Several other physical examples are worked out in 
Appendix O 

To be able to add a Majorana mass at the y = boundary, we need of course to 
consider a fermion that belongs to a real representation of the unbroken gauge group. On 
top of the bulk action ()2.1|) . we then consider the following brane action (slightly separated 
from the boundary): 

S 4D = J d A x\L(M X X + M* xx) ly=e • (4-7) 

Note that the mass has been written as ML to give to M a mass dimension equal to one. 

To find the modified BC's in the presence of this brane mass term, the first step is 
to chose the Dirichlet-Neumann BC's the two two-component spinors \ and ip would have 
satisfied at the "real boundary" y = in the absence of the Majorana mass term. Thus as 
previously we assume that: 

<9 5 X|o = 0, V|o = 0, d 5X \ L = 0, if>\ L = 0. (4.8) 

The effect of the brane mass term is to introduce discontinuities in the wave functions at 
y = e. The bulk equations of motion are modified to 

-ia^d^x ~ d 5 ip + rrMfj + .1/7. 6{y - e) x = 0, (4.9) 
-%o*d$ + d 5X + mx = 0. (4.10) 

Integrating the first equation over the delta function term shows that, while x remains 
continuous, the value of the ip profile undergoes a jump: 

% e = M*Lx\e- (4.11) 
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Because ip undergoes a jump, the second bulk equation of motion requires that the deriva- 
tive of x also undergoes a jump: 



to] 



= lO' 



(4.12) 



In the limit e — ► and from the fixed values (J4.8|) of \ an d ip at y = 0, the jump equations 
finally give the BC at y — + . This is what we will be interpreting as the BC corresponding 
to the theory with Majorana masses on the boundary y = (the BC at y = L remain of 
course unaffected by the mass term localized at y = 0): 

<9 5 X|o+ = i^d^\ + - mX\o+, = MLx\o+, d 5 x\L- = 0, 4>\ L - = 0. (4.13) 

The first equation is just the bulk equation of motion evaluated at the boundary, so we 
conclude that the BC corresponding to the Majorana mass is just 



What one should recognize at this point is that (|4.1Hjl corresponds precisely to the boundary 
condition 



with ML = —sq/cq. Thus we have a completely dynamical description of one of the 
boundary conditions mentioned in Section El One can also reproduce the other types of 
BC's by adding different boundary localized operators like a localized kinetic term or a 
mass interaction term with boundary localized fermions. An exhaustive list of cases are 
worked out in Appendix O 



5 Fermion Masses in the SU(2) L xSU(2) i? xU(l) B _ L Flat 
Space Toy Model 



As an application of the previous sections, we will consider generating masses for leptons in 
a higgsless extra dimensional model in flat space where the bulk gauge group ^ is SU(2)/,x 
SU(2)^x U(1)_b_l as presented in Ref. pQ. This model has large corrections to the elec- 
troweak precision observables, and can not be viewed as a realistic model for electroweak 
symmetry breaking. However, most of the large corrections can be eliminated by putting 
the same model into warped space 0, or as pointed out recently in [TSj by adding large 
brane localized gauge kinetic terms in the flat space case. We find it useful to first present 
some of the features of the construction for the fermion masses in the flat space model, as 
a preparation for the more complicated warped case presented at the end of this paper. 

• With quarks in the bulk there is also a bulk SU (3) co ior gauge group. 



lp\ + = MLx\o+- 



(4.14) 



(coip + s x)\o+ = 0. 



(4.15) 
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5.1 Lepton sector 



As always in a left-right symmetric model, the fermions are in the representations (2, 1, —1/2) 
and (1,2,-1/2) of SU(2) L x SU(2) fi x \J(1) B _ L for left and right handed leptons respec- 
tively. Since we assume that the fermions live in the bulk, both of these are Dirac fermions, 
thus every chiral SM fermion is doubled (and the right handed neutrino is added similarly). 
Thus the left handed doublet can be written as 

(Xu L ,$u L ,Xe L ,$e L Y , (5.1) 

where (Xu L ,Xe L ) will eventually correspond to the SM SU(2) L doublet and (ip UL ,ip eL ) is 
its SU(2) l antidoublet partner needed to form a complete 5D Dirac spinor. Similarly, the 
content of the right handed doublet is 

(Xu R ,^u R ,Xe R ,^e R Y , (5.2) 

where (ip VR) ip eR ) would correspond to the 'SM' right-handed doublet, i.e., the right electron 
and the extra right neutrino, while (Xu R ,Xe R ) is its antidoublet partner again needed to 
form a complete 5D Dirac spinor. 

For simplicity we assume for now the absence of a bulk Dirac mass term (we will later 
need to introduce such terms in the warped scenario to build a realistic model). In the 
absence of any brane induced mass terms the fields Xu L , Xe L , ip VR and ip eR would contain zero 
modes, while the other fields would acquire a KK mass of the order of the compactification 
scale. Thus without the additional boundary masses that we need to add, the boundary 
conditions would be 

^L\0,L = ^eL\0,L = Xvr\o,l = Xe fl |o,L = ®- (5-3) 

We also assume that, as explained in [BE!, on one brane the SU(2)lX SU(2)# symmetry 
is broken to the diagonal SU(2)d, while on the other brane SU(2)#x XJ(1)b-l is broken to 
U(l)y. This means that on the SU(2)d brane the theory is non-chiral, and a Dirac mass 
term Md connecting the left and right fermions can be added. Assuming that this Dirac 
mass term is added on the y = L brane, the boundary conditions will be the same as in 
Appendix IC.2I 

*Pe Llo+ = 0, X e Rl o+=Q, (5-4) 
^e L \ L - = -M D Llp eRlL _, Xe R \ L - = M D L X e L \ L -- (5.5) 

The KK decomposition is of the form (|3.1|) . leading to an electron mass being solution of 
the equation: 

tan(m n L) = M D L (5.6) 

which, for <C l/£, is solved by m ~ M^. The lowest mass state is as expected a Dirac 
fermion with a mass that is just given by the Dirac mass added on the brane. 
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However, the unbroken SU(2)^ symmetry so far guarantees that the neutrino has the 
same mass as the electron. The neutrino mass needs to be suppressed by some sort of a 
see-saw mechanism, which can be achieved, as in Appendix IC.3[ by coupling the neutrino to 
a fermion localized on the brane where SU(2)^x U(1)b-l is broken to U(l)y. Let us thus 
introduce an extra right-handed neutrino £ VR localized on that brane. Being SU(2) L xU(l)y 
neutral, this extra brane fermion can have a Majorana mass as well as a mixing mass term 
with ij) VR via the 4D Lagrangian at y = 

5 4D = J d A x ( - itu^d^vn + ML 1 ' 2 (i v J) VR + £uJ> VR ) + f (UU + CuJu R )j • ( 5 - 7 ) 

The boundary conditions on the SU(2)lX U(l)y brane are then 

^ L , 0+ =0, (5.8) 
X u R \^ = -ML^% R , (5.9) 
-i^d^XvR^ + fXu R \o+ - M 2 L ipv R \o+ =0. (5.10) 
The boundary conditions on the SU(2)d brane remain untouched: 

W| L - = -M D Li> VRlL _, (5.11) 
Xu R \ L -=M D L X u L \ L -. (5.12) 
The KK expansion will be of the form 

£ = $>„£„(£), (5-13) 

n 

Xv L = Y.9n L \z)U*), Xv r = Y,9 { : r \*)Ux), (5-14) 

n n 

= Yl fn L \ Z ) Znix), = E fn R Hz) UX), (5-15) 

n n 

where the £ n 's are 4D Majorana spinors of mass m n : 

-ia^d^n + m n £ n = and - ia^d^ n + m* £ n = 0. (5.16) 

Together with the bulk equations of motion these BC's lead to the following mass spectrum 
for neutrino's (assuming the eigenmass m n real) 

2 M2L n {MpL cos 2 m n L - sin 2 m„L) = (1 + M 2 D L 2 ) sin2m„L. (5.17) 

For f'MpL <C M 2 and M^L < 1 we get that the lowest mode is a Majorana fermion with 
a mass approximately given by 

mo ~ (5.18) 

which is of the typical see-saw type since the Dirac mass, Mp, which is of the same order as 
that of the electron mass, is suppressed by the large masses of the right handed neutrinos 
localized on the brane. Thus a realistic spectrum is achievable in this simple toy model for 
the leptons. 
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5.2 Quark sector 



In order to get a realistic mass spectrum for the quarks, one cannot simply add a single 
brane localized two-component fermion as for the neutrinos, since in that case we would 
induce an anomaly in the effective theory. Instead, we need to introduce a vector-like brane 
localized color triplet with the quantum numbers of the up-type right handed quark and 
its conjugate (or the down type for a mixing for the down quarks). So the fields that we 
are considering now are 



/ Xu L \ 

Xd L 



Xur ^ 
Xd R 



Vur 



(5.19) 



\ i>d L J \ ^d R ) 



where (Xu L ,Xd L ) will be identified as the SM SU{2) L quark doublet and tp UR and ipd R as 
the SM right handed quarks; (ip UL , ipd L ) and Xu R and Xd R are their partners needed to 
form complete 5D spinors and they will get KK masses of order of the compactification 
scale. Finally (£, UR ,Vu R ) is a localized 4D Dirac spinor that will couple to tp UR on the 
SU(2)l x U(l)y brane at y — 0. We will thus again assume that in the absence of the 
brane localized mass terms and mixings the fields \u L , Xd L , ^u R and ipd R would have zero 
modes. The 4D brane localized terms at y = are: 



y=0 



d 4 x 



+ / {VuJur + VuJur) + ML 1 ' 2 (£ UR <lf> UR + CuJu R )^ , (5.20) 

while at y = L, we just had an SU(2)^ invariant Dirac mass term mixing the left and the 
right quarks: 



S 



d*x M D L {Xu L ipu R + Xu R ^u L + h.c.) + (Xd L ipd R + Xd R ipd L + h.c.) (5.21) 



y=L 



Since we have not included any mixing term at y = for the down-type quarks their 
spectrum will just be of the same form as for the electrons above, determined by the 
equation 

tan(m„L) = M D L. (5.22) 

The boundary conditions for the up-type quarks are similar to those obtained in Ap- 
pendix IC.31 



VVl|o+ — 0' 
X« fl |o+ = —ML 1 / 2 £ UR , 
(d^ + f 2 )Xu Rl0+ ~ M 2 Lia»d^ URl0+ = 0, 
^ L \l- = ~M D L tp UR |£ _ , Xu R \ L - = M D L Xu L \ L - 



(5.23) 
(5.24) 
(5.25) 
(5.26) 
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The KK decomposition will be of the usual form ()3.1|) and leads to the following quantiza- 
tion equation similar to the neutrino's mass equation: 

mp — f 2 

2 -f- 2 — — (sin 2 m n L - M 2 D L 2 cos 2 m n L) = (l + M 2 L 2 ) sin 2m n L, (5.27) 

which again, for / M^L <C a// 2 + M 2 and M^L <C 1, has a solution approximated by 

Mpf 

m ~ — :. (5.28) 

If M -C /, we can suppress the brane Dirac mass using a see-saw type mechanism, or if 
M ~ / we can get a mass just slightly modified compared to the brane Dirac mass. We can 
use this freedom to generate both the masses of the light first two generations as well as 
the masses of the massive third generation. We will discuss this much more in the context 
of the more realistic warped model in the coming section. 



6 Fermions in Warped Space 
6.1 Bulk equations of motion 

We now extend the discussion of the previous section to a truncated warped spacetime, 
which as usual we take to be a slice of AdSs [8] . The conformally flat metric corresponding 
to this situation is given by 

ds 2 = (^pj (rj^dx^dx" - dz 2 ) . (6.1) 

The boundaries of the spacetime are at R ~ 1/Mp/ and R' ~ 1 TeV -1 . Fermions in such 
a space have been considered in jTUHnillZllIHllini ■ Here we first briefly review the generic 
features of the fermion wave functions in this space, and then repeat our analysis for the 
acceptable boundary conditions for this situation. 

The fermion action in a curved background is generically given by 

S = J d 5 x^ ef ^DmV - D M V ef 7 a tf ) - m^j , (6.2) 

where is the generalization of the vierbein to higher dimensions ( "fiinfbein" ) satisfying 

e£W = 9 MN , (6-3) 

the 7 a 's are the usual Dirac matrices, and Dm is the covariant derivative including the spin 
connection term. Note again that the differential operators have not been integrated by 
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parts in order to avoid the introduction of any boundary terms that would otherwise be 
needed to make the action real. 

For the AdSs metric in the conformal coordinates written above, e a M = (R/z)51 I , and 
= (d/j, + 7^75/ (4^))^, = d$ty, however the spin connection terms involved in the 
two covariant derivatives of ()6.2j) cancel each other and thus do not contribute in total to 
the action. Finally, in terms of two component spinors, the action is given by 



,S* = / <l '.r I - 



v ix^d^x-i^d^ + l(^x-X^^) + ^(^X + X^ , (6-4) 

where c is the bulk Dirac mass in units of the AdS curvature 1/R, and again d§ = d§ — d§ 
with the convention that the differential operators act only on the spinors and not on the 
metric factors. 

The bulk equations of motion derived from this action are 

- ia^x ~ + = 0, (6.5) 

z 

-io*d$ + d 5X + "—^X = 0. (6.6) 

The KK decomposition takes its usual form (the case of a Majorana KK decomposition 
will be discussed in details in Section 16.2. 2j) 



X 



5^ 9n(z)Xn(x) and ip = ^ f n (z)4> n (x), (6.7) 



where the 4D spinors Xn and ip n satisfy the usual 4D Dirac equation with mass m n : 

-ia^d^Xn + m n ili n = and - ia^d^ n + m n x n = 0. (6.8) 

The bulk equations then become ordinary (coupled) differential equations of first order for 
the wavefunctions /„ and g n : 

fL+mn9n-^fn = 0, (6.9) 

c-2 

g n -m n g n A — g n = 0. (6.10) 

For a zero mode, if the boundary conditions were to allow its presence, these bulk 
equations are already decoupled and are thus easy to solve, leading to: 

c+2 



fa = Co (-) , (6.11) 



z 



2-c 



9o = A , (6.12) 
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where A$ and Co are two normalization constants of mass dimension 1/2. 

For the massive modes, the first order differential equations can be uncoupled and we 
obtain the two second order differential equations: 

fn ~ ifn + ( m l ~ f ^ ± )fn = 0, (6.13) 
9'n ~ % + i< ~ C ^)9n = 0, (6.14) 

whose solutions are linear combinations of Bessel functions: 

g n (z) = z'i (A n J c+ i(m n z) + B n Y c+ i(m n zj) (6.15) 
f n (z) = z* (c n J c _^{m Tl z) + D n Y c _i(m n z)^j . (6.16) 

The bulk equations of motion ()6.9|) - (j6.10J) further impose that 

A n = C n and B n = D n . (6.17) 



6.2 Boundary conditions 

To find the consistent boundary conditions, we need to again consider the boundary terms 
in the variation of the action: 

R' 

(6.18) 

R 

which agrees with the expression for flat space up to the irrelevant factor of R 4 /z 4 . Thus the 
boundary conditions that make the boundary variation of the action vanish will generically 
be of the same form as for the flat space case, that is those given in equations ()2.8|) - (J2.11J) . 
If the fermions are in real representations of the gauge group M£ is allowed to be non- 
vanishing, but if they are in complex representations M@ has to be zero. 



fiSb oun d — — I d X 



R 



6.2.1 Boundary conditions in absence of extra boundary operators 

As an example let us consider the simplest case, when we make the conventional choice of 
imposing Dirichlet BC's on both ends: 

i> lR+ =0 and i) ]RI - =0. (6.19) 

These BC's allow for a chiral zero mode in the x sector while the profile for ip has to be 
vanishing, so we find for an arbitrary value of the bulk mass c that the zero modes are 
given by: 

/ = and go = A ^y~ C . (6.20) 
16 



The main impact c has on the zero mode is where it is localized, close to the Planck brane 
(around z = R) or the TeV brane (around z = R'). This can be seen by considering the 
normalization of the fermion wave functions. To obtain a canonically normalized 4D kinetic 
term for the zero mode, one needs 



f - ) o A l ( o ) = 1 *- e - = - ^== (6.21J 



z J R \RJ RWR n - 2c -R 



where the first factor in the integral comes from the volume ^/g, the z/R factor from the 
vierbein and the rest is the wave function itself (squared). To figure out where this zero 
mode is localized in a covariant way, we can send either brane to infinity and see whether 
the zero mode remains normalizable. For instance, sending the TeV brane to infinity, 
R' — > oo, the integral ()6.21j) converges only for c > 1/2, in which case the zero mode is 
localized near the Planck brane. Conversely, for c < 1/2, when the Planck brane is sent to 
infinity, R — > 0, the integral (j6.21|) remains convergent and the zero mode is thus localized 
near the TeV brane. In the AdS/CFT language ^2] this corresponds to the fact that for 
c > 1/2 the fermions will be elementary (since they are localized on the Planck brane), 
while for c < 1/2 they are to be considered as composite bound states of the CFT modes 
(since they are peaked on the TeV brane). 

This result can also be seen easily when using the proper distance coordinate along the 
extra dimension. In this case the AdSs metric is written as (k = 1/R is the AdS curvature): 

ds 2 = e~ 2ky rj^dx^dx" - dy 2 . (6.22) 
And the actual normalized wavefunctions (including the volume and vierbein factors) are 

e -(2c-l)k{y-y Pl ) for c > y 2 ( and _^ ^ 

e -(l-2c)k(y TeV -y) for c < ( and _^ _ oc ^ (Q 2A j 

Finally, let us point out that if we were to impose Dirichlet BC's on both ends of the 
interval for x, we would have found a zero mode in the ip sector. And this zero mode would 
have been localized on the Planck brane for c < —1/2 and localized on the TeV brane for 
c> -1/2. 



6.2.2 Boundary conditions with a brane Majorana mass term 

To familiarize ourselves more with the BC's in warped space, we will repeat the flat-case 
analysis of section [Ol and consider the case when a Majorana mass is added on the Planck 
brane for the x field (which would otherwise have a zero mode).-'' Based on our discussions 
we expect that the boundary condition on the Planck brane would be modified to 

(cosa^ — sin a x)| R + = (6.25) 

^While this paper was in preparation |19| appeared, which also presents a detailed treatment of a 
Planck-brane localized Majorana mass term. 
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where sin a = corresponds to the case with no Majorana mass, while cos a = to the case 
with a very large Majorana mass. To identify the actual relation between the Majorana 
mass and a we again consider adding the Majorana mass at z = R + e, read off the BC's 
from the bulk equations, and then send e — > + . In this case the bulk equation of motion 
will be modified to 

c + 2 - M*R 2 

-ia^d^x - + ——tp + 6(z-R-e) = 0, (6.26) 

where M is the Majorana mass added. There will be a discontinuity in the profile for if), it 
is given by the jump equation 

l4>] lR+e = MR X \ R+ e, (6.27) 
from which, using ip\ R = 0, we can read off the relevant boundary condition 

i)\ R+ =MR X \r + , (6.28) 

which is indeed of the form ([6.25)1 . 

The KK decomposition is of the form 

X = ^2gn(z)£ n (x) and = ^ fn(z)£ n (x), (6.29) 

n n 

where the 4D spinors £ n satisfy the usual 4D Majorana equation with mass m n : 

-ia^d^n + m n l n = and - ia^d^ n + m*£ n = 0. (6.30) 

Instead of the expansion in the Bessel and Neumann functions, it turns out that it is more 
convenient to expand in terms of J v and J^ v . These functions are linearly independent as 
long as v is not an integer, that is if c 7^ 1/2+ integer. We will then treat the c = 1/2 as a 
special case later. The reason why it is more convenient to use these functions is that the 
expansion for small arguments, needed for an approximate solution for the lowest modes, 
will be much simpler if we use this basis. Thus the wave functions /„ and g n will be of the 
form 

g n (z) = z'i (A n J c+ i(\m n \z) + B n J_ c _i(\m n \zj) , (6.31) 
f n (z) = z% (c n J c _i (\m n \z) + D n J_ c+ i(\m n \z)^j . (6.32) 

and the bulk equations further require that 

m n A n = \m n \C n and m n B n = -\m n \D n . (6.33) 
The two boundary conditions 

^ + = MR X \r + and = (6.34) 
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then lead to the equation determining the eigenvalues m n : 

J c -J- c+ i 2 - J-c+Jc-i = ±MR{J c+ iJ_ c+ i + J_ c _i J c _i), (6.35) 



where J u = J u (\m n \R) and J u = J v (\m n \R'). This equation can be approximately solved 
for the lowest eigenmode (assuming that m R' <C 1) by expanding the Bessel functions for 
small arguments as 

(TT tv 1 7v ( 6 - 36 ) 



.i) r(z/ + i)' 

We find, that for c > 1/2 the lowest eigenmode is approximately given by 

m ~(2c-l)M, (6.37) 

while for c < 1/2 it is 



m ~ (1 - 2c)M ^-j . (6.38) 

The c = 1/2 case has to be treated separately, since in that case the expansion has to be 
in terms of the Bessel and Neumann functions. The equation that the eigenvalues have to 
solve will be given by 

J Y -Y J = ±MR(J 1 Y -Y 1 J ). (6.39) 
For the lightest mode we find 

M 

m o ~ r^r- 6.40 

The interpretation of these expressions is quite clear. When c > 1/2, we expect the resulting 
mass to be proportional to the mass added on the Planck brane, since the fields themselves 
are localized near the Planck brane. For c < 1/2 the zero mode is localized near the 
TeV brane, so adding a mass on the Planck brane has only a small effect due to the wave 
function suppression. For the c = 1/2 case the wave function is flat, and one expects the 
usual volume suppression as in flat backgrounds; that is one expects a suppression by the 
proper distance between the branes. The expressions above are in clear correspondence 
with these expectations. 



7 Fermion Masses in the Higgsless Model of Elec- 
troweak Symmetry Breaking in Warped Space 

We are now finally ready to consider the SU(2)ixSU(2)ijxU(l)B-L model in warped space, 
where electroweak symmetry breaking is achieved by boundary conditions (rather than by 
a Higgs on the TeV brane). As discussed in 0, this model has a custodial SU(2) symmetry 
that protects the p parameter from large corrections, and thus to leading log order the 
structure of the standard model in the gauge sector is reproduced. However, an obvious 
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lingering question is whether realistic values of the fermion mass could be obtained in this 
model. This model could be considered the AdS dual of walking technicolor and it is 
well-known that in technicolor theories it is difficult to naturally obtain a realistic fermion 
spectrum. Here we will show, that in the extra dimensional model there is enough freedom 
in the parameter space of the theory to be able to incorporate the observed fermion masses. 

7.1 Quark sector 

The left handed leptons and quarks will be in SU(2)x doublets, while the right handed 
ones in SU(2) fi doublets, exactly as in (j5.1|) . f)5.2|) and (|5.19|) . Thus we will have two SU(2) 
doublet Dirac fermions for the leptons and two separately for the quarks in the bulk for 
every generation, xl,r, $l,r- Each Dirac fermion has a bulk mass cl,r and a Dirac mass 
Md on the TeV brane that mixes the two bulk fermions. In addition we assume that there 
is a Dirac fermion localized on the Planck brane that mixes with ipR. This will again be 
necessary to be able to sufficiently split the masses of the up and down type fermions. 

We again assume that the fields i/)l,r and Xl,r are such that in the absence of brane 
localized masses/mixings the fields xl and i/)r would have zero modes, that is the BC's in 
the absence of the brane terms are as in ()5.3j) 

^l\r ]R > = Xr\r,r> = 0. (7.1) 

Since we would like the zero modes (at least for the light fermions) to be localized near 
the Planck brane (in order to recover the SM relations for the gauge couplings), we need 
to pick the bulk mass terms cl > 1/2 and cr < —1/2. The reversal of the inequality for 
cr is due to the fact that for the right handed doublets we want the ip fields to have zero 
modes, and for these types of zero modes the localization properties as a function of c are 
modified, with c < —1/2 localized near the Planck brane while c > —1/2 near the TeV 
brane. 

The bulk part of the fermion action will be as in ()6.2|) . and the bulk equations of motion 
are as in (j6.5|) - ()6.6j) . To find the appropriate boundary conditions, we need to consider the 
brane localized mass and mixing terms. The mixing term on the Planck brane will be of 
the form 

S Pl = Jd 4 x(^- iZ<r»d4 - i V <T"dtf + f(vC + + MVR + &r^J ^ , (7.2) 

where £ and rj are brane localized fermions, which together form a Dirac fermion with a 
Dirac mass / on the brane. On the Planck brane only SU(2)^xU(l)y is unbroken, and 
this extra Dirac fermion is assumed to be an SU(2)^ singlet carrying the U(l)y quantum 
numbers of the right handed SM fermion fields, such that the mixing with mixing mass 
M in (I7.2|) is allowed. This is the analog of the y = term in the Lagrangian in the flat 
space case discussed in Sections 4 and 5. Since the warp factor on the Planck brane is one, 
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XR\R+ = m n— J2—?pR\R+, (7-5) 



the boundary conditions following from this brane localized Lagrangian will exactly match 
that in flat space: 

(d^ + f)xR lR+ = M 2 Ria^ RlR+ , (7.3) 
ipL lR+ = 0. (7.4) 

For modes with m\ <C f 2 the BC ()7.3|) can be approximated by 

M 2 R 
1 

where m n is the n th mass eigenvalue. 

The mass term on the TeV brane will be given by 

S T eV = J M D R' (^ rXl + Xl^r + i>LXR + Xr^lJ • (7.6) 

Note, that in order to have a natural Lagrangian the parameter M^R' should be of order 
one, thus should be of the order of TeV. From this mass term the boundary conditions 
on the TeV brane will be analogous to the flat space case discussed in Sections 4 and 5: 

iI>l\r?- = -M d R!iJ> R \ B ,- (7.7) 
XR\Rf- = M D Bf X L\Bf- (7-8) 

For the mode decomposition the bulk wavefunction solutions we take for the general 
case 1/2 + cl )R ^ integer: 

Xl,r = A\Al R Jx {m n z) + Bl R Jx ( 

*Pl,r = A[A n LR J_i (m n z)-Bl R Ji (m n z)), (7.9) 

V ' 2^ L ' R 2 L ' R J 

where m n is the 4D mass of the given mode that one is considering. For a mode with 
m n R' < 1 we can expand the Bessel functions for small arguments, and since the coefficients 
AL,R,Bi tR depend on the eigenvalue m n , some overall powers of m n can be absorbed into 
these constants to make the expansion more transparent (from here on we will suppress 
the index n), keeping the terms at most quadratic in m n , we get: 

. A LR mz CL ' R+1 ~ nr ( m 2 z 2 

XL,R = Z 2 — + Br R Z- C ^ 1 



2c L)R + l ' l V 2-4c LjR 



*lr = z 2 [A LR z^[l- mV )_ B ^mz^ R 
1 LM 1 2 + 4c Lii? ; l-2c L;R 1 
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Imposing the above boundary conditions we find that the lightest eigenmode is approxi- 
mately given by§ (assuming again cl > 1/2, c R < —1/2) 



y/(2e L -l)(-2c R -l)fM D /iA"- 1 



V/ 2 + (-2c fl - 1)M 2 U 

The approximation cl > 1/2, cr < —1/2 should be sufficient for the light quarks, but for the 
top quark we need the wavefunctions to have a larger overlap on the TeV brane. Therefore 
we need to consider the case cl < 1/2, c R > —1/2 so that the top quark is localized near 
the TeV brane rather than near the Planck brane, which will make it possible to get a 
sufficiently large top quark mass. Using similar methods as before we find that in this case 
( cl < 1/2, c R > —1/2) the approximate lowest mass eigenvalue is 

m 2 M 2 (1-2c l )(1 + 2cr) 

° 1 + MIR*(1 -c L + c R ) + (§ ) 2c * +1 Ml (i + 2CR ) (l + " 



For completeness we also briefly discuss the special cases when cl = 1/2 and/or 
c R = 1/2. With cl = 1/2, c R > —1/2 we find that the lightest eigenmode is given 
by 

m2 2{l + 2c R )pMlR>^ 

m ° ~ log f ( J R /1 + 2c «/ 2 (2 + (1 + 2c R )M 2 D R' 2 ) + 2,R 1 + 2c *((l + 2c H )M 2 - / 2 )) 1 ' ' 

The complementary case c R = —1/2 and cl > 1/2 gives 



(2c L -l)/ 2 M 2 fR\ 2cL - 1 
/ 2 logf + M 2 UV 



-o - \,:^ R /, ^ [^) > (7-14) 



while for the doubly special case cl = —c R = 1/2: 



2 / 2 ^D / x 

m ° ~ logf^logf + M 2 ) ' (7 - 15) 



Let us now use these expressions to demonstrate that it is possible to obtain a realistic 
mass spectrum for all the SM fermions. We will use (J7.11|) for the first two generations, 
while for the third generation quarks we use (|7.12j) . We have also numerically solved the 
bulk equations with the appropriate boundary conditions, and found that ([7. 11)1 - ([7. 12)1 are 
generically good approximations for the lowest eigenvalues, up to the ten percent level. 
For the results presented below we have used the numerical solutions to the eigenvalue 
equations rather than the approximate formulae (J7.11j) - (J7.12)) . We also use the numerical 

§In the following formulae we only keep the terms which can be leading in R 1 / R for chosen values of 
bulk masses. Among the remaining terms we separately keep only the leading contributions proportional 
to M 2 and f 2 terms because M/f is a free parameter which may be large. 



22 



solution to find the lightest KK excitations in each case. We will not attempt to explain all 
the observed CKM matrix elements in this paper, though we see no reason why it should be 
hard to obtain the right values. In order to correctly reproduce jS] the masses of the W and 
the Z, throughout the fits we will use the values R = 1(T 19 GeV" 1 , R' = 2 • KT 3 GeV -1 . 
The parameter Mo should be of the order of the TeV scale, while the splitting between the 
up and down-type fermions will be obtained by choosing an appropriate value for the ratio 
Mj f, that is the ratio of the mixing with the brane fermions to the diagonal mass of the 
brane fermions. A reasonable quark mass spectrum can be obtained using the following 
parameters: 

(i) For the first generation take cl = —cr = .6, Mo = 50 GeV and M/f = 3.8 for the 
up sector and M/f = for the down sector. Then m u w 3 MeV, « 6 MeV and the 
first KK excitations appear at 1.2 TeV then 1.3 TeV (both for up and down). 

(ii) For the second generation take cl = — cr = .52, Mo = 112 GeV and M/f = 50 
for the strange sector. Then m s ~ 110 MeV, m c w 1.3 GeV and the first KK excitations 
appear at 1.1 TeV then 1.3 TeV (both for s and c). 

(iii) For the third generation we need localize both the left handed and the right handed 
zero modes near the TeV brane in order to be able to get a large enough top quark mass. 
One numerical example is c L = 0A,c R = -1/3, M D = 900 GeV, / = 2.5 ■ 10 10 GeV and 
M = 10 15 GeV for the bottom sector and M = for the top sector. For these parameters 
we get mtop ~ 175 GeV and my, pa 4.5 GeV. The first KK excitations of the bottom quark 
appear at the relatively low value of ~ 550 GeV, while for the top quark at ~ 700 GeV. 
This would imply that the third generation (since it would be localized near the TeV brane) 
would be very different from the first two, and interesting effects in flavor physics could 
be observable. For a recent analysis of examples of the consequences for a composite third 
generation see |25] . 

The above numbers are only given for the purpose of demonstrating the viability of 
obtaining a realistic set of quark masses. However, there are several free parameters that 
one can vary for obtaining the correct masses: cl,cr and Mo, while the ratio M/f is 
mostly set by the amount of splitting within a multiplet. Here we have only assumed the 
simplest possibility when one of the two fermions within a generation have a mixing with the 
brane localized fermions. Clearly, there is a much richer spectrum of possibilities for such 
mixings which we will however not deal with here. To illustrate some of the available free 
parameters, we have varied the c's and Mo around the solution for the second generation, 
while keeping the c and s quark masses fixed. The resulting relation between cl,cr and Mo 
is displayed in Fig. Ufa). In Fig. C^b) we show the dependence of the mass of the lightest 
KK mode on the parameters cl and cr. Note, that a characteristic feature of all of these 
solutions is that the mass of the lightest KK mode decreases with increasing Mo- This is 
due to the fact that for large Mo the KK mass is sensitive to the value of the c's but not 
Mo itself. This can be simply seen by taking the large Mo limit, where the resonance mass 
is just set by the scale 1/R'. 

Since the amount of interesting new flavor physics in the third generation crucially 
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Figure 1: (a) Contour plot of the value of needed to obtain m c = 1.2 GeV forvarying 
values of Cl and Cr. The regions starting with the darkest moving toward the lighter 
ones correspond to M D = 0.1,0.3,0.6,1,1.5,2,2.5,3,3.5,4,4.5 TeV. (b) Contour plot of 
the value of the lightest KK mass for the second generation quarks assuming that Md is 
chosen such that m c = 1.2 GeV, for varying values of cl and cr. The regions starting with 
the darkest moving toward the lighter ones correspond to Mkk = 0.1,0.3,0.5,0.7,0.9, 1.1 
TeV 

depends on the deviation of from 1/2, we have examined how small cl — 1/2 could be. 
For this we have generated calculated the acceptable values of cl,cr and Md that would 
give us the correct top quark mass, which is shown in Fig.EJ One can see, that the smallest 
possible value for 1/2 — cl is ~0.03. 

It is also interesting to note that since the mixing with very heavy fermions on the 
Planck brane is essentially equivalent to introducing brane kinetic terms, we can easily 
implement the Hiller-Schmaltz mechanism for solving the strong CP problem If all 
intergeneration mixing arises on the Planck brane by mixing with the heavy Planck brane 
fermions, then the net effect for the light fermions is that all the mixing appears in kinetic 
terms, then all complex phases can be rotated into the CKM matrix without introducing 
strong CP violation [2l)] . 

7.2 Lepton sector 

One has a variety of options for generating the lepton masses. The nicest possibility would 
be to have an extra dimensional implementation of the usual neutrino see-saw mechanism 
which takes advantage of the fact that the right handed neutrino is a singlet under all 
the SM gauge groups. This implies that on the Planck brane one can simply add a brane 
localized mass term to the right handed neutrino (since SU(2)# is broken there). Using 
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Figure 2: Contour plot of the value of Mo needed to obtain m t0 p = 175 GeV for varying 
values of cl and cr. The regions starting with the darkest moving toward the lighter ones 
correspond to Mb = 1, 1.5, 2, . . . , 6, 6.5 TeV. 



similar methods as before we find that the neutrino mass is given by 

Ml / R\ 2 ^-"- 1 ' 



where Mr is the Majorana mass of the right handed neutrino on the Planck brane. In 
absence of any brane localized fermions, the charged lepton masses are given by (|7.1H) 
with M = then we find the relation between the neutrino and the charged lepton mass 
to be 

1 m r 

m Vi = — (7.17) 

-2c R -1 Mr 

where the m,- are the masses of the charged leptons. Note, that this is almost completely 

i 

analogous to the usual see-saw formula, with the only difference being that the large scale 
directly suppresses the charged lepton mass squares, and not a mass of order 100 GeV as is 
usually assumed. That is, the difference is in the appearance of the charged lepton Yukawa 
coupling. Using this mechanism one could get realistic lepton masses using for example the 
following parameters: 

(i) For the first generation take c L = -c R = .65, M D = 130 GeV and M R = 10 10 GeV. 
Then m e « 500 keV, m Ue m 10~ 7 eV. 

(ii) For the second generation take Cl = —cr = .58, Mp = 250 GeV and Mr = 
3 ■ 10 10 GeV. Then m M w 100 MeV, w 2 ■ 10~ 3 eV. 

(iii) For the third generation take cl = —cr = .53, M D = 240 GeV and Mr = 
10 12 GeV. Then m T w 1.7 GeV, m Vr w 4 ■ 10" 2 eV. 

Instead of the appearance of the intermediate scale Mr ~ 10 11 GeV for the Majorana 
mass of the right-handed neutrino (which is somewhat lower than usually assumed), one 
can insist on the real see-saw formula. A real see-saw mechanism can actually be achieved if 
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Figure 3: (a) Contour plot of the value of Mr needed to obtain m M = 100 MeV and 
m v = 2 • 10 -3 eV for varying values of Md and M/f. The regions starting with the darkest 
moving toward the lighter ones correspond to Mr = 4, 6, 8, 10, 13, 16, 19, 22 • 10 14 GeV. (b) 
Contour plot of the mass of the first KK excitation of the muon neutrino, keeping fixed m M = 
100 MeV and m v = 2 ■ 10~ 3 eV and varying Md and M/f. The regions starting with the 
darkest moving toward the lighter ones correspond to m* K = 300, 400, 500, 600, 700 GeV. 

the suppression of the charged lepton masses compared to Md is the consequence not only 
of a warp factor suppression but also of an additional mixing with Planck-brane localized 
fermions as in the general case (|7.11|1 . The neutrino masses would still be given by ([7.16)1 . 
with a scale Mr different from the scale M in (|7.11|) . This way one could choose an Md 
of order .1 to 2 TeV for all three generations, and the scale Mr needed for the neutrino 
masses would be closer to the usually assumed values of order 10 15 to 10 16 GeV. A possible 
choice of parameters could be for example: 

(i) For the first generation take c L = -c R = .55, M D = 100 GeV, M/f = 1500 for the 
suppression of the electron mass to m e pa 500 keV, and Mr = 10 16 GeV for the suppression 
of the electron neutrino mass to m Ve ~ 10~ 8 eV, 

(ii) For the second generation take c L = -c R = .52, M D = 1000 GeV, M/f = 500 
for the suppression of the muon mass to pa 100 MeV, and Mr = 10 15 GeV for the 
suppression of the muon neutrino mass to m v pa 2 ■ 10 -3 eV, 

(iii) For the third generation take c L = -c R = .51, M D = 2000 GeV, M/f = 100 
for the suppression of the tau mass to m T pa 1.7 GeV, and Mr = 6 • 10 14 GeV for the 
suppression of the muon neutrino mass to m Vr pa 3 ■ 10 -2 eV. 

On Fig. Efa) we have plotted the values of the see-saw scale, Mr, needed to reproduce 
= 100 MeV and m v = 2 ■ 10~ 3 eV, while varying the other free parameters, Md and 
M/f (we have further imposed that cl = —cr). On Fig. Efb), we have also computed 
the the mass of the first KK excitation of the muon neutrinos. Its mass decreases as Md 
increases and we note that it is also almost independent of the value of M/ f needed to fit 
the muon as long as we keep the mass of the muon neutrino fixed. 
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Another possibility that we will not explore here, due to its relative complexity, is that 
the charged lepton and neutrino come from two different SU(2)# doublets as in ref 

8 Conclusions 

We have considered theories with fermions on an extra dimensional interval. We derived the 
consistent BC's from the variational principle and explained how to associate the various 
BC's to different physical situations. We have applied our results to higgsless models of 
electroweak symmetry breaking, and showed that realistic fermion mass spectra can be 
generated without the presence of a Higgs fields both in flat and in warped space. 
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Appendix 

A Spinor and Gamma Matrix Conventions 

For completeness, we give in this appendix the convention about spinors and Dirac ma- 
trices used throughout the paper. We have mainly followed the conventions of Wess and 
Bagger [20J. 

We are working with a mostly "— " space-time signature, (H ), and we have 

chosen the chiral representation of the Dirac gamma matrices: 



( 



a" 




) 



/i = 0,1,2,3 and T 5 



( 



il 2 






-il 2 



) 



(A.l) 
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where and <r M are the usual Pauli matrices 

'-(:;)•••-(! oWro. (A2) 

A famous relation about the Pauli matrices is 

a»a v + a v a» = 2rT and a»a u + a u a» = 2^ . (A.4) 

A 5D Dirac spinor is written in terms of a pair of two component spinors 

* = I $ I (AJS) 



The dotted and undotted indices of a two component spinor are raised and lower with the 

2 

i/3 



2x2 antisymmetric tensors e a p = i/r^p and = ia^- and their inverse e a/3 = —icr^p, 



a/3 

X a = e^XfS and $ & = e &$ ^. (A.6) 

Note also the adjoint relation: 

(X f ) a = t- (A.7) 
Finally, x>P and X^ denote the two Lorentz invariant scalars: 

X^ = X a ^ a and x4> = Xa^° '■ (A.8) 
These products are symmetric 

X^ = ^X and X^ = ^X- (A.9) 

B Examples of KK Decomposition in Flat Space with 
the Simplest BC's 

We present in this appendix explicit examples of KK decomposition of fermions in flat 
space. 

We begin by giving the full KK decomposition in the case of the simplest Dirichlet- 
Neumann BC's: there are a priori four different cases to discuss, but the case X|o — X\l = 
is similar to ip\ — i>\L — and the case x\o — = is similar to ip\o — X\l = 0. So let us 
mention the KK decompositions for the latter cases only: 
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(?) when ip\ = ip\i = 0: there is a zero mode for x oiu Yj with an exponential wave 
function localized either on the or L brane depending on the sign of bulk Dirac mass, as 
well as a tower of massive modes mixed between \ and ip: 

oo 

X = A e~ my £ + ^ A n (cos k n y - —sin k n y) Xn, (B.l) 



, k n 

n=l 



00 



ip = -^2~^A n smk n y ijj n , (B.2) 

n=l 11 

with KK masses m n solution of the quantization equation: 

sin k n L = 0, k 2 n = m 2 n - m 2 . (B.3) 

(ii) when x/j\ = x\l — 0: there is no zero mode at all and the tower of massive modes 
is given by 



X = ^ A n (cos k n y - —sin k n y) Xn, (B.4) 



CO 

m 



ip = -^2—^A n sink n y ip n , (B.5) 

n=l n 

with KK masses m n now solution of the equation: 

k n = m tan k n L, k 2 n = m 1 n — m 1 . (B.6) 

Let us now discuss the KK decomposition when we impose BC's of the form ()2.12|) : 

So,l i/j\o,l + c ,l X\o,l = (B.7) 

where so,l (co,l) stand for the sine (cosine) of some (possibly complex) angles, ao,L, that 
determine which linear combination of the fields on the two boundaries are vanishing, the 
first thing to notice is that these BC's can be imposed only when the 5D fermion belongs to 
a real representation of the gauge group. In this case, the KK modes will be 4D Majorana 
fermions and the KK decomposition will be of the form 

X = Y,9niv) Ux), (B.8) 

n 

^ = X)/n(3/)£n(aO. (B.9) 

n 

with the spinors £ n { x ) satisfying the Majorana equation: 

-^«9^ + m n £ n = 0, (B.10) 
-i^d^n + ml^^O. (B.ll) 
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The wave functions are of the form (|3.9|) - (|3.1U|) with m n being now replaced by \m n \. In 
general there is no zero mode except if there exists a tuning between the angles cko.l and 
the bulk mass m: 

sin(a;o — a£) — sin(ao + (Xl) tanhmL = 0. (B.12) 

For the massive KK modes, as before, the bulk equations give A n and B n in terms of C n 
and D n and the boundary equations reduce to two complex equations for two complex 
unknowns. The masses are thus obtained as the roots of a 4 by 4 determinant. After some 
algebra we obtain the quantization equation: 

(|soSl| 2 + |coCl| 2 ) |m 2 1 sin 2 k n L — \sQCi,\ 2 {rn sin k n L — k n cos k n L) 2 
-\c s L \ 2 (m sin k n L + k n cos k n L) 2 + {c s* s L c* L + s c* c L s* L ) k 2 n = 0. (B.13) 

When all the boundary angles are real, the above equation factorizes into the two simpler 
equations: 

±m n cos(«o — oil) sin k n L — k n sin(a — oll) cos k n L + m sin(a + «l) sin k n L = 0. (B.14) 
In particular, we recover the results ()B.3|) - (jB.6|) for the particular angles a , «l = 0, tt/2. 



C Boundary Conditions in the Presence of Various 
Localized Mass/Mixing/Kinetic Terms 

In this appendix, we would like to extend the discussion of Section HJ and present the BC's 
in the presence of various localized operators on the boundaries. Of particular interest 
are localized terms for fermions in complex representations. Majorana masses for such 
theories are forbidden, however, localized kinetic terms are not, neither are localized Dirac 
masses. We will also discuss examples where bulk fermions are mixing with fermions that 
are localized at the boundaries. 

We will be following the steps outlined at the beginning of Section HJ we first add the 
localized terms at y = e to the bulk equations of motion in order to efficiently deal with 
the discontinuities in the wave functions. At the real boundary y = we impose as always 
the simplest BC's 

<9 5 x|o = 0, V>|o = 0, <9 5X |l = 0, ip\ L = 0. (C.l) 
We then take e — > and identify the relevant BC's at y = + and y = L~ . 



C.l Adding localized kinetic terms 

We add an extra boundary kinetic interaction localized at y = e > to one of the two 
two-component spinors through the 4D boundary action 

5*40 = - / d^xmx^d^Wy^- (C.2) 
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The equations of motion are modified to 

-ia^x (1 + K6(y - e)) - + mtfj = 0, (C.3) 
-ia^d^ + d 5 x + mx = 0. (C.4) 

Integrating the first equation over the delta function gives a jump condition for ip which 
then implies a jump in the derivative of x- 

[ip] |e = -ma^d^e and [<9 5 x] k = [$] |e . (C.5) 

Using the values (|C.1|) of the fields on the boundary, we finally obtain the BC's we were 
after 

^| + = -ma M <9 At X|o+ and <9 5 X|o+ = i^d^^ . (C.6) 

This is exactly a boundary condition of the form (-0 — ica^d^x) |o+ = obtained from the 
variational principle in (|2.11j) . where the parameter c being identified as the coefficient of 
the boundary localized kinetic term. 



C.2 Adding a Dirac mixing of bulk fermions on the boundary 

We consider two bulk fermions that are mixed together through a Dirac mass on the 
boundary 

S^d = j d A x M D L ^VaXi + Xi4>2 + ^1X2 + X 2 ^iJ (C.7) 

We further need to specify the values of the fields at y = and y = L. We are assuming 
that as usual 

d 5Xl\0,L = °> ^l|0,i = °> d 5^2\ ,L = 0) X2| ,L = 0, (C.8) 

in such a way that when Mr> — > there are two zero modes corresponding to Xi and fa- 

In the present case, there is an ambiguity when we want to push the interaction mass at 
a distance e away from the boundary Indeed if we were to push the whole expression (|C.7j) . 
we will end up with discontinuity in the wave functions of both ipi and fa and Xi and X2, 
which then requires a regularization of the products like fa5(y ~ e), for instance by an 
averaging of fa over its limits from both sides. In this case, using the values (1C.8J) of the 
fields at y — 0, we get the following BC's at y = + 

/ MdL 

^ + = l + iM^ V ' ( } 

M D L 

X2 ' 0+ -~1 + ±M^ Xl ' 0+ - (C - 10) 
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Another regularization would consist for instance in smoothing the delta function by a 
square potential, M D /aQ(y)Q(a — y), and then take the limit a — ► 0. In this case, we 
would arrive at slightly modified BC's of the form 

V>i|o+ = 2tanhM D L/2 ^ 2 | +, (Cll) 
X2| + = -2tanhM D L/2 Xi|o+- (C.12) 

Another possibility will be to impose the values (jC.8|) on the boundary before pushing 
it e away from the brane. In that case, the bulk equations of motion will be compatible 
with continuous wave functions for ip 2 and Xi an d would require discontinuities in ipi and 
Xi only. No further regularization of the delta function would then be necessary. And the 
BC's would simply read 

^i| 0+ =M D Lij 2l0+ , (C.13) 
X 2\o+ = ~M D L X i\o+- (C14) 

We note that whichever way we go, the BC's are all the same in the limit of small 
MoL. And, in the flat case, the lowest eigenmode is a Dirac fermion with mass Mq. More 
importantly, the form of the boundary condition is independent of the regularization of the 
delta function, and at most the interpretation of the physical meaning of Mr> might depend 
on it. 



C.3 Mixing with brane fermions 

We now analyze a more general case where, in addition to the bulk fermions, there are 
localized spinors on the boundaries which mix with bulk fermions through Dirac mass 
terms. The addition of a Dirac mixing mass on the y = L brane leads to a 4D action given 
by 

S iD = Jd 4 x(^- i&»d4 - i W ^fj + ML 1 ' 2 (nx + XV) + f + £fj) J ^ , (C15) 

where M and / have mass dimension 1. 

To find the corresponding boundary conditions, we follow our by now usual procedure 
and push the interactions a small distance e away from the boundary at y = L, and solve 
the resulting equations of motion which are given by 

-ia^d^x ~ + m^ + ML 1/2 fj 5(y - L + e) = 0, (C. 16) 

-ia"d li $ + d B x + mx = 0, (C17) 

-w»dfi + ML l ' 2 x + /£ = 0, (C. 18) 

-i^d^ + fTj = 0. (C.19) 
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Integrating the first equation over the delta function gives the jump for ip: 



= ML 1 / 2 fj, (C.20) 
which, using the defined value (|C.1J) of tj) at y = L implies the following BC at y = L~ 

^ |L _ = -ML 1/2 fj. (C.21) 
Then the last two equations can then be combined in the following way: 

+ f) fj = -iML 1 ' 2 a^ X \L. (C.22) 
Taking the limit e — > 0, we finally arrive at the BC's: 

= -ML^fj, (C.23) 

(d^ + f 2 )^ = iM 2 La^ X \L- ■ (C.24) 



We notice that interactions (jC.15j) lead to a generalization of the boundary condition 
obtained by adding a brane localized kinetic term. Indeed, in the limit when / and M tend 
to infinity with a fixed ratio, the boundary condition (jC.24|) reduces precisely to ()C.6|) . This 
is expected since in the limit of large mass / of the brane fermions they must be integrated 
out. This leaves one massless fermion, which is a linear combination of £ and \i an d its 
kinetic term contains a delta-function contribution, originating from the localized kinetic 
term for £, which is proportional to M 2 /f 2 . When we constructed a realistic pattern for 
the quark and lepton masses in Section [?1 we explicitly introduced brane localized fermions 
but it should be kept in mind that a similar spectrum can be obtained just by introducing 
on the Planck brane localized kinetic terms for the bulk fermions. 

When the bulk fermion appears to be neutral with respect to the residual gauge sym- 
metry on the boundary, it can couple to a Majorana fermion on the brane. So let us see 
how the previous BC's are modified in that case. The 4D boundary localized action is 

S 4D = [ d 4 x( - irta^fj + ML 1 ' 2 (r] X + XV) + lf(w + Vv)) (C.25) 

J V " J \y=L 

where M and / have again a mass dimension 1. 

The same procedure of pushing the interaction e away from the boundary and taking 
the limit e — > leads to a jump in ip and a jump in d$x an d the BC's finally read: 

fy L - = -ML 1/2 fj, (C.26) 
d,x\L- = {-M 2 L X + m\L-. (C.27) 
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